Analysis of measurement errors for a superconducting phase qubit 
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We analyze several mechanisms leading to errors in a course of measurement of a superconducting 
flux-biased phase qubit. Insufficiently long measurement pulse may lead to nonadiabatic transitions 
between qubit states |1) and |0), before tunneling through a reduced barrier is supposed to dis- 
tinguish the qubit states. Finite (though large) ratio of tunneling rates for these states leads to 
incomplete discrimination between |1) and |0). Insufficiently fast energy relaxation after the tun- 
neling of state |1) may cause the repopulation of the quantum well in which only the state |0) is 
supposed to remain. We analyze these types of measurement errors using analytical approaches as 
well as numerical solution of the time-dependent Schrodinger equation. 

PACS numbers: 85.25.Cp, 03.67.Lx, 74.50.+r 



I. INTRODUCTION 

In recent years, a significant progress has been made 
in developing superconducting Josephson-junction cir- 
cuits for quantum computation. Single-qubiliiiS* and 
coupled-qubi1i^i& operations have been demonstrated 
for various qubit types, including charge?^ flux^ and 
phas o 2 ! 5 ! 6 ! 7 ! 8 ! 9 qubits. Recently introduced flux-biased 
phase qubit a 6 i 7 i 8 i 9 (Fig. 1) significantly reduce the total 
dissipation during measurement in comparison with the 
current-biased phase qubits^ which leads to lower de- 
coherence and faster recovery after measurement. 

Measurement of qubits is an important stage in quan- 
tum information processing. A novel scheme for fast 
measurements of flux-biased phase qubits has been 
demonstrated— and subsequently employed for simulta- 
neous measurement of two coupled qubits^ According to 
this scheme [see Fig. 1(b)], a measurement is performed 
by lowering the barrier between the qubit "left" potential 
well and a much deeper "right" well, so that tunneling 
from qubit state |1) (the first excited left-well state) to 
the right well occurs with probability close to one, while 
state |0) (the left- well ground state) remains intact. 

It is of interest to analyze the measurement process 
theoretically, since such an analysis can be used to min- 
imize measurement errors. In measurements of coupled 
qubits, there are generally one- and two-qubit errors. The 
main source of two-qubit errors is crosstalk^ arising af- 
ter the application of the measurement pulse due to the 
coupling between the qubits. The crosstalk was analyzed 
by us recently^ 

The present paper is devoted to the theoretical study 
of the behavior of a single flux-biased phase qubit during 
the measurement pulse, which lowers the barrier. Our 
purpose is to find potential sources of one-qubit mea- 
surement errors in the fast measurement scheme used in 
Refs. and 0, and analyze relations between the error 
level and the parameters of the measurement pulse. We 
consider three mechanisms of single-qubit errors. First, 
there is a possibility of transitions between qubit states 
|1) and |0) due to nonadiabaticity of the process. Second, 
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FIG. 1: (a) Circuit schematic of the flux-biased phase qubit, 
controlled by external magnetic flux 4>{t). (b) Sketch of the 
qubit potential energy U(S), where 5 is the superconducting 
phase difference across the Josephson junction. During the 
measurement pulse the energy barrier is lowered so that the 
state |1) escapes from the "left" well into the neighboring 
"right" well; this escape is sensed by a nearby SQUID. 



there may be an incomplete discrimination between the 
two qubit states by tunneling. Third, it happens that if 
the after-tunneling energy relaxation is not fast enough, 
there is a possibility of a repopulation of the left qubit 
well (in which only the state |0) is supposed to remain) 
after tunneling of the state |1) to the right well. 

The paper is organized as follows. Section ITU contains 
a discussion of flux-biased phase qubits, measurement 
errors, and the method for numerical solution of the 
Schrodinger equation used in the paper. Nonadiabatic 
effects are studied in Sec. IIIII where we derive analytical 
formulas and perform numerical calculations. The effi- 
ciency of the qubit state discrimination by tunneling is 
discussed in Sec. IV. The repopulation error mechanism 
is discussed in Sec. V, where we perform full numerical 
simulation of the quantum evolution in absence of energy 
dissipation, and estimate the rate of dissipation required 
to eliminate the repopulation effect. Finally, Sec. VI 
contains concluding remarks. 



II. MODEL 

The flux-biased phase qubit [Fig. 1(a)] has the same 
circuit schematic as the rf SQUID iii In its quantum de- 
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scription the wavefunction *f?(5,t) evolves according to 
the Schrodinger equation with the Hamiltonian 



2m 



U(5,t), 



(1) 



where S is the superconducting phase difference across 
the Josephson junction, p = —ih(d/dS) is the momentum 
operator, m = [$o/(27r)] 2 C is the effective mass, C is the 
junction capacitance, <£>o = h/(2e) is the flux quantum, 
e is the electron charge, and the qubit potential is 



U(6,t)=Ej 



2A 



cos 5 



(2) 



where Ej = I^q/{2tt) is the Josephson energy, A = 
27r/ L/$ is the dimensionless inductance, <j)(t) — 
27r<i>(t)/$o is the dimensionless external magnetic flux, 
&(t) is the magnetic flux in the loop, In is the critical 
current, and L is the inductance. We will assume that 
before the measurement pulse the right well of the qubit 
potential [s ee F ig. 1(b)] is much deeper than the left well 
(as in Refs.|a7-8.9). The qubit levels |0) and |1) are, re- 
spectively, the ground and the first excited levels in the 
left well. 

The depth AUi >r of a left (I) or right (r) well (i.e. the 
energy difference between the potential maximum and 
minimum) depends on external flux <fi and can be char- 
acterized by the estimated number of the discrete levels 
in the well: 



AUi 



l.r 



Hull r 



(3) 



where uj^ r is the classical oscillation frequency near 
the well bottom (the "plasma frequency"), uj^ r = 
y/Ej{l/\ + cos 8i tr )/m (here 5i :T corresponds to the well 
bottom). The energy dissipation in the phase qubit can 
be described by introducing resistance R into the circuit 
[Fig.^a)]; however, for the most of this paper the dissi- 
pation is neglected (R = oo). 

The fast qubit state measurement is achieved by ap- 
plying a pulse of the magnetic flux <j>, which lowers suffi- 
ciently slowly (adiabatically) the barrier height AUi (and 
so Ni), until the first excited state |1) of the qubit left 
well becomes very close to the barrier top, so that the 
system in this state can easily tunnel through the bar- 
rier to the right well, while tunneling from the ground 
state |0) is almost negligible. (In principle, the measure- 
ment can be arranged in a way so that the state |1) goes 
over the barrier; however, in this case tunneling from the 
state |0) becomes significant.) After the measurement 
pulse the potential U(S) is returned to its initial shape 
by decreasing <fi back to its value before the measurement. 
Ideally, the system initially in the upper qubit state |1) 
should switch after the measurement to the right well, 
while the qubit initially in the ground state |0) should 
remain in the left well. However, at least two kinds of 
errors are possible in this process of qubit measurement. 
First, finite duration of the measurement pulse (which is 



supposed to be rather fast) and corresponding nonadia- 
batic effects during increase of the flux 4>(t), may lead to 
transitions between levels |1) and |0) (and other levels) 
before the tunneling starts. Second, during the tunneling 
stage the qubit in the state |1) may erroneously remain 
in the left well and/or state |0) may erroneously switch 
to the right well. One reason for the second type of error 
is finite ratio of tunneling rates for the states |1) and |0) 
(ideally, it should be infinitely large). Another reason is 
the incomplete switching of the state |1) due to repopula- 
tion of the left well from the right well during or after the 
tunneling stage. These kinds of errors will be considered 
in the following sections. 

To solve numerically the time-dependent Schrodinger 
equation, we use the multi-projection approach of Ref. 
ha In more detail, we divide the pulse duration r 
into small time increments At and at each moment 
tk = kAt (k = 0, 1, . . . ) we find the eigenvalues E* and 
eigenf unctions ipn(S) for the Hamiltonian H(tk), using 
the Fourier grid Hamiltonian method^ (equivalently, the 
periodic pseudospectraU^ method) . Then the wavefunc- 
tion is computed as 



t) = £ a*^(5) e- 1 ^-^ 1 (t k <t< t k+1 ) 

n 

where 



(4) 
(5) 



is calculated using the wavefunction ^(S, tk) at the end 
of the previous timestep. Starting the evolution from the 
initial qubit level i, we calculate the evolving population 
of the level n as 



Pni{t) 



k |2 



for t k <t < tk+i- 



(6) 



III. NONADIABATIC EFFECTS 



During the stage of rising flux <f> (before the tunnel- 
ing), the population of the state |0) is supposed to re- 
main unchanged. (Notice that excitation of state |1) to 
higher levels does not lead to the measurement error.) 
However, since increase of <f> is not infinitely slow, there 
may be nonadiabatic transitions between levels |0) and 
|1) (and/or higher levels) leading to measurement errors. 
Such nonadiabatic errors are discussed in this section. 

For numerical calculations we use the circuit parame- 
ters of Ref. la 



C = 700 fF, L = 0.72 nH, In = 1.7 mkA. 



(7) 



Also, we assume that during the first half of the measure- 
ment pulse the left-well dimensionless barrier N (here 
and below we omit the subscript I) decreases from N = 5 
to = 1.355, which corresponds to <fi varying from 5.09 
to 5.31 (for TV = 1.355 the state |1) is very close to the 
barrier topi^). 
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FIG. 2: The modified qubit potential (SJ (thick solid lines, in 
GHz) and the normalized wavefunctions |^ n (<$)| 2 (thin solid 
lines), shifted vertically by the energy eigenvalues E„ for (a) 
N — 5 and (b) N — 1.355; other qubit parameters are given 
by J7J. Dashed lines correspond to the actual potential 
The energy origin is chosen at the well minimum. 



For simplicity, in this section we neglect tunneling. 
Correspondingly, we consider a modified potential energy 
(see Fig. , which differs from Eq. (J2J) by absence of the 
right well, so that the potential is constant to the right 
of the barrier-top position S b : 



U(6,t) 



Ej{(6- 
U(6 b ,t), 



) 2 /2A-cos6}, 5<5 b , 
S > 6 b , 



(8) 



where both <fi and 8 b change with time. In numerical 
calculations we limit the range of S between 0.5 and 3, 
so that the wavefunction ^>(5,t) is assumed to vanish at 
S = 0.5 and 6 = 3 (we have checked that further increase 
of this range does not change the results). Figure[3shows 
the wavefunctions |"0n(^)| 2 and energies E n for the states 
localized in the well, corresponding to the beginning and 
the maximum of the measurement pulse. Since the nona- 
diabatic transitions after the tunneling stage are not im- 
portant, in this section we consider only the first part of 
the measurement pulse (increase of 4>, decrease of N ) . 



A. Analytical theory 

1. Nonadiabatic error 

Assuming that the qubit wavefunction before the mea- 
surement is 



|*o> = c \ip ) +ci\ipx), 



(9) 



we describe the qubit evolution during the rising part 
of the measurement pulse by expanding its wavefunction 
over the instantaneous basis given by the cigenfunctions 
ijj n (S,t) of H(t) with the eigenvalues E n (t): 

= Y,Cn{t)e^l n ^o E ^ dt 'M5,t)- (10) 



and 



the initial conditions are co(0) = Co, ci(0) 
c„(0) = for n > 2. 

Let us define the nonadiabatic measurement error Q 
as the difference between the ground state occupation 
Po(t) — l(V'o(i)| v I'(i))| 2 before the measurement pulse 
(t = 0) and at the maximum of the pulse (t = r): 



Q = |Po(0)-P o (r)|. 



(11) 



Here Po(0) — |co| 2 , while Pq{t) can be expressed via the 
evolution operator S(t) [so that = S(t)|*o)]: 



Po(t) = |co| 2 Poo(t) + | Ci | 2 P 01 (t) 
+2Re[cocJ5oo(r)5 * 1 (T)], 



(12) 



where S ni (r) = (V>„(t)|S(t)|^(0)), while P m (r) = 
I'S'ni ("J") | 2 is the occupation of state n for initial state i. 
The last term in the rhs of Eq. (|12fl results from quantum 
interference. 

We are interested in the case of a small nonadiabatic 
error, Q <C 1. Then, as will be shown below, Sqi(t) and 
1 — 5*00 (t) are quantities of the first and second order in 
the perturbation, respectively. Hence, up to the second 
order, the nonadiabatic error is given by 

Q = I M 2 [l - Poo(r)] - |ci| 2 Poi(r) - 2Re[c c* 1 S* 1 (T)}\, 

(13) 

Notice that 1 — P 00 and Poi are both of the second order 
in perturbation; therefore in the case when both qubit 
states are initially occupied, the first two terms in Eq. 
(|13H can be neglected in comparison with the third (in- 
terference) term. The third term is maximized when both 
qubit states are equally occupied, |co| = |ci| = \/\/2, and 
coc*5'q 1 (t) is a real number. So, the maximum nonadia- 
batic error is approximately 



\S i (r) | - y/^Jj). 



(14) 



Notice that if initial state is l^i), then the nonadiabatic 
error is Q = P i(r). If initial state is \ipo), then the error 
is 1 — Poo(t) = Xm>i Pno( T ), which as will be seen later 
is approximately equal to Poi(t) (because Pro ~ Poi in 
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the lowest order, while P n o for n > 2 are very small 
because direct transitions between non-neighboring levels 
are suppressed). Therefore, if only one qubit state is 
initially occupied, the error (in the lowest order) is 



Pc = Pqi(t) = QI 



(15) 



which is much smaller than Q e . This fact can be easily 
understood by visualizing operator S as a rotation of 
the qubit Bloch sphere (neglecting higher levels). Then 
rotation by a small angle (p (around £-axis) leads to the 
change of z-coordinate by ip near equator (for x — z — 0) , 
while the change is only ip 2 /2 near the poles. 

In this section we will mainly analyze the error P e ; 
however, it should be remembered that the actual nona- 
diabatic error can be up to Q e = \fP~ e . 



2. Adiabatic perturbation theory. 

The evolution operator S(t) obeys the Schrodinger 
equation, which can be reduced to the following 



(xp n (t)\H{t)\^ m {t)) cl ^ L 

TiUnm. {t} 



^ dt 'S mi) (16) 



where u> nm — (E n — E m )/h. (The state label i appears 
only in indices, while i in the exponent is the imaginary 
unit.) The initial condition for Eq. {Tfj} is 5„j(0) = 5^, 
where S^ is the Kronecker symbol. 

Equation (|16l) is exact. It can be simplified assuming 
slow variation of H(t), so that H can be treated as a small 
perturbation. In the zero approximation (the adiabatic 
approximationiSiii) the rhs of Eq. 1|16H is assumed to 
vanish, so that S — and therefore S n i(t) = 8^. In the 
first approximation we use S m i = S^ in the rhs of Eq. 
H16| l. which yields^ 7 (for n ^ i) 



hu> ni (t') 



(17) 

Notice that for n 7^ i, S n i(t) is of the first order in per- 
turbation, while 1 — Sa(t) vanishes in the first order [see 
Eq. Ijlfip]. so that non-zero 1 — Su(t) appears only in the 
second order, when the first order S m i is used in Eq. (|16|) . 
Also notice that in the first order S n i(t) = —S* n (t), as 
follows from Eq. (|T7|l ; therefore P ni (t) = P in (t) . We have 
used these properties in the previous subsection. 
For the modified potential JSJ we find 



H(t) = U(S,t) 



Ej[4,(t)-S]<p(t)/X, 5<6 b (t); 
dU(S b (t),t)/dt, 5>6 b (t). 



For simplicity let us use the formula in the first line of Eq. 
(I18H even at S > 8 b (this approximation can be justified 
by the fact that the wavefunctions are small to the right 



of the barrier top position S b )- Then the matrix elements 
of H can be calculated as 



(il> n (t)\H(t)\TPi{t)) - -Ejj>{t)6 ni {t)/\, (19) 



where 



(20) 



is the "position" matrix element. Inserting Eq. pity 
into l|17fl . we obtain the following equation for P n i(t) = 
\Sni{t)\ 2 (n 7^ i) in the lowest order of perturbation: 



Pni{t) 



El 



, <f>{t')5 ni {t') i( t> U)ni(i » )di ' 



(///— — — ' 



.( 21 ) 

Notice that for purely harmonic potential the matrix el- 
ements Sni are non-zero only for levels which are nearest 
neighbors. This means that P n i for non-nearest neigh- 
bors would appear only in higher orders in perturbation. 
Since the qubit potential is not purely harmonic, there 
will be direct transitions between non-neighbors; how- 
ever, these transition are significantly weaker than for 
nearest neighbors. 

When the measurement pulse duration is scaled with- 
out change of the pulse shape, one can write the external 
flux as 



0(t) -> (b(t/r) 



(22) 



where r is the duration of the rising part of the pulse and 
6 — t/r. Then the Hamiltonian and hence its eigenstates 
and eigenvalues depend only on 8: 

H(t)^H(9), LJ ni (t) -> <j ni (6), iMt)-»V»(0). (23) 
As a result, Eq. I|21|) can be recast in the form 



Pni{t) 



El 



del t'(0')Sni(e') ciTf ? Uni( e'> )de » 

[i UJni{0') 



(24) 

where </>'(6>) = dcj)(9)/d6. 

Finally, the nonadiabatic measurement error defined 
as P e — Pqi{t) [see Eq. lfT5| ] is obtained as 



El 
hX 



uw(0) 



(25) 



Even though this is only the lowest order approximation 
for Poi(r), the numerical calculations of Pol show (see 
below) that Eq. (|25|) is quite accurate. 

This formula for P e can be further simplified (though 
with a noticeable loss of accuracy) by considering the 
qubit as a harmonic oscillator with the frequency aj w 
changing in time. In this approximation^ 



(26) 



5 



so Eq. I|25(l becomes 



E 2 j 



2h\ 2 i 



dd 



iT/ fl Wl o(e')de' 



(27) 



This expression requires only the knowledge of the in- 
stantaneous transition frequency ujioit/r) an d does not 
depend on the instantaneous wavefunctions. 

Even further simplification is possible if we neglect the 
change of the frequency lo\q with time (as we see later, 
this is a quite crude approximation). Then the mea- 
surement error P e is simply determined by the spectral 
component of the derivative 0' of the external flux at the 
frequency ujiq: 



Pr 



E] 



2h\ 2 



<t>'{9) e iuJl0T9 d6 



(28) 



Since u>io changes during the measurement pulse quite 
significantly (by more than 50 %, see Fig.[2Jl, application 
of Eq. I|28[) is not straightforward and requires a choice of 
the time moment within the pulse, at which wio is taken. 
However, Eq. I|28|) can still be used for simple estimates of 
the nonadiabatic error, for example using the maximum 
and minimum values of wio (at the beginning of the pulse 
and at the end of its rising part). Notice that Eq. i|28|) 
implies that the nonadiabatic error can be significantly 
reduced by suppressing the spectral components of (j)(t) 
in the frequency range of wio- 



3. Long-r behavior. 

In this subsection we analyze the asymptotic behaviour 
of the nonadiabatic measurement error P e for long pulse 
duration r, assuming that the pulse profile <f)(t) is an- 
alytic everywhere except for the endpoints t = and 
t = t. We show that in this case the long-r behavior 
of P e can be expressed via pulse properties at the sin- 
gular (nonanalytic) endpointsiS Even though the result 
obtained below in this subsection cannot be directly ap- 
plied to a realistic experiment, it still sheds a light on 
the dependence of the measurement error on the pulse 
shape. 

The j-fold integration by parts of Eq. ("P7fl yields 



S n i{r) 



e Jo 



+i J 



i [p u ni (t')dt' d 3 

dP 



dt k - 



k-l 



(*) 



e J o 



g(t) 



dt, 



o 

(29) 



where g(t) = (tp n (t)\H (t)\ipi(t)} /~hiv n i(t) and the level la- 
bel i again appears only in indices (i in other places is the 
imaginary unit). For scaling the pulse duration r without 
change of a pulse shape, we use relation d/dt = T~ 1 d/d9 
[see Eqs. ("2"3) and """j"}]; then the fcth term in Eq. (j!"9"|) 



acquires a factor r~ k . Let us denote by jo and j\ the 
lowest orders of the time derivatives of the Hamiltonian, 
which do not vanish at t = and r, respectively. Then 
for sufficiently large r, (mincj„i)r 3> 1, one can leave in 
Eq. H29|) only the lowest-order nonvanishing terms cor- 
responding to t = and r (9—0 and 1) and obtain the 
approximation 



S n i(T) 



i*>(^ n (o)|ff<*>>(o)|^(o)> 



K5 +1 (0)rA 
i^ n (l)|ffO-0(i)|^(i)) 



TlLJ 



ji+i 



(l)rii 



(30) 



where time is normalized by r, H^' (9) is the jth deriva- 
tive of H{9), and uj n i — J io n i{9)d9 is the frequency 
averaged over the pulse duration. 

Using the relation (yj n {9)\H^ {9)^(9)) 
-Ej(f>^(9)d m (9)/X, which is similar to Eq. fKty . 
we obtain the following approximation for the nonadia- 
batic error P e = |Soi(t)| 2 : 



E 2 j 



h 2 \ 2 T 2 i 



(jo) (O)<5io(O) 



Jo- 



_1 (0)tJo- 



iii-jo0ai)(i)5 lo (i) 



ji+i 



(1)t* 



(31) 



where j = minjjo, ji}. The error probability l|31l) gen- 
erally decays with r as t~ 2j . In particular, in the case 
jo = ji = j, Eq. 1(3*11) becomes 



P, 



E 2 j 



h 2 \ 2 T 2 i 



0k'>(O)<MO) 0W(l)*io(l) 



(0) 



(1) 



(32) 

In this case besides decaying with r, the dependence 
P e (r) also oscillates with the frequency ljiq. In contrast, 
when jo 7^ ji > on e of the two terms in Eq. (|31() dominates 
at large r, and the decay occurs monotonously. 



4- Time- dependent perturbation theory 



The approach used in subsection IIII A 21 can be com- 
pared with the standard time-dependent perturbation 
theory. Let us split the Hamiltonian (TJ into the "un- 
perturbed" Hamiltonian Hq and "perturbation" V(t): 



H(t) = H 
U (S) = E, 



V(t), Ho = 

(6 - <t> r ? 
2A 



P 

2m 



cos S 



V{t) = -{E J /\)[<j ) {t)-(j> r 



(33) 

(34) 
(35) 



where we neglect the change of the energy origin 
and choose an arbitrary flux <fi r as a reference point. 
Even though the standard time-dependent perturbation 



6 



theory^ can be only applied when the perturbation V(t) 
is small, which is obviously not the case in our problem, 
let us still apply it formally. Then for the nonadiabatic 
error we obtain 



p = f Ej5 w 
\h\ujxo 



= ( EjSw 
\h\ujxo 



(j)(t)e luJl0t dt 



(j)'{6)e iuJl0Te d9 



(36) 



where <5io and luiq correspond to the Hamiltonian H . It 
is easy to see that Eq. I|36|) coincides with Eq. (|25|l if the 
time-dependence of S\q and wio in Eq. (|25[1 is neglected. 
An arbitrary choice of <p r in the unperturbed Hamiltonian 
with a natural limitation 4>(0) < <fi r < (/>(r) corresponds 
to an arbitrary choice of the time moment within the 
pulse, at which S\o and wio are taken. 

Equation (|36|l can be simplified by using the approx- 
imation (|26p. Then it becomes exactly Eq. i|28[l . which 
expresses P e via the spectral component of the derivative 
of the external flux at the transition frequency loiq. Even 
though this simplest result can be obtained in both ap- 
proaches, we would like to emphasize that the approach 
of subsection IIII A 21 allows us to use much more accurate 
approximations for the error probability P e . 



B. Numerical results 

For numerical calculations we use the method de- 
scribed in Section II [see Eqs. which corresponds 
to the solution of exact equation Q16[l. Let us describe 
the time dependence of the flux during the rising part of 
the measurement pulse as 



<K*/r) = (f> + {(t>x - <h) X g(t/T), 



(37) 



where 4>o and aH are the initial and final flux values, while 
g(9) describes the pulse shape and satisfies conditions 
g(0) = and g(l) = 1. Unless mention otherwise, we 
choose in this subsection <fio — 5.09 (corresponding to 
N = 5) and cj)\ = 5.31 (corresponding to N = 1.355). 

Consider first the simplest case of a linearly increasing 
pulse, g(9) = 9. Figure [3] shows the time evolution of 
three quantities: Poii -PlOj an d 1 — -Poo (which all could 
in principle be used for the definition of the nonadiabatic 
error) for two values of the pulse duration: r — 0.2 ns 
and r = 2 ns. One can see that as expected, these three 
quantities practically coincide: 



P i(*)«-Pio(i)«l--Poo(t), 



(38) 



so that it is not really important which of them is used for 
the definition of the error P e . (In this subsection for cal- 
culation of P e we will use 1 — Poo> which is slightly larger 
than two other quantities.) The dominating process is 
transitions between states |0) and |1); as a result, the 
oscillations in Fig. [3| are at the frequency lo\q (t) [see Eq. 
(|2ip]. The oscillation amplitude increases with time in 
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FIG. 3: Time dependence of the instantaneous nonadiabatic 
error defined in three different ways: Poi(i) (solid lines), 
-Pio(i) (dashed lines), and 1 — Poo(i) (dotted lines), for two 
durations of the linear pulse: r = 0.2 ns and r = 2 ns (for 
r = 2 ns the data are multiplied by the factor 20) . The solid, 
dashed, and dotted lines are practically indistinguishable, ex- 
cept near the pulse end for r = 0.2 ns. 



Fig. |21 This is explained by the fact that increasing flux 
during the measurement pulse decreases the transition 
frequency a>io (see Fig. |5J , which leads to the increase of 
the integrand in Eq. lt2*T|l for Poi(t) [note also Eq. (|2l)|) ]: 
therefore, the oscillation amplitude should increase ap- 
proximately as [wio(i)] -3 - Figure |3| also demonstrates 
that the instantaneous nonadiabatic error is much larger 
for t = 0.2 ns than for r = 2 ns. 

Figure 0] shows the dependence of the error P e on the 
duration r of the linear pulse on the double-logarithmic 
scale. The nonadiabatic error generally decreases with 
increase of r (for linear pulse crudely as t -2 ); however, 
the dependence also shows oscillations, which originate 
from oscillations in Fig. |3| In Fig. 01a) the thick solid 
line shows the numerical result. It is practically indis- 
tinguishable from the thin solid line showing the lowest- 
order analytical result (|2~5ll . The dashed line corresponds 
to the simplified formula (|27|l : one can see that it is also 
quite close to the exact result. We have observed a similar 
relation between the results obtained numerically, by an- 
alytical formula l|25|l . and by simplified formula l|27|l for 
other pulse shapes discussed below. Figure 0^b) shows 
the dependence of the error on the pulse duration for 
three values of the final barrier height: N = 1.355, 1.5 
and 1.65 (the initial value is N — 5 for all the cases). One 
can see that the error increases with the decrease of the 
final barrier height; however, dependence is not strong 
and therefore the results presented here are not quite 
sensitive to the exact choice of the pulse amplitude. 

Figure [S] shows the comparison between the numeri- 
cal results for P e (r) dependence and the results obtained 
using Eq. I|28|l which relates P e to the spectral compo- 
nent of 4> at the frequency luiq. Three panels of Fig. 
are for (a) the linear pulse, (b) the pulse with shape 
g{9) = 1 — cos 4 (7r6'/2), and (c) Gaussian pulse. (For 
Gaussian pulse we define r as the r.m.s. width of the 
pulse; the fact that the pulse is really longer than r, 
reduces P e in comparison with other shapes.) The fre- 
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FIG. 4: The nonadiabatic error P e as a function of the dura- 
tion r of the linear measurement pulse, (a) Numerical results 
(thick solid line) and results obtained using Eq. I|2~>t (thin 
solid line) and Eq. 1271 (dashed line) for the pulse starting 
with the barrier height N = 5 (corresponding to <f>o = 5.09) 
and ending with N = 1.355 (01 = 5.308). (b) Numeri- 
cal results for pulses starting with N = 5 and ending with 
N = 1.355 (thick solid line), 1.5 (dashed line), and 1.65 (thin 
solid line); corresponding maximum fluxes are <j> = 5.308, 
5.298, and 5.287. 
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quency wio in Eq. (|28|l is not well-defined since it changes 
during the pulse by more than 50%; so in Fig. we show 
two analytical results, corresponding to at the be- 
ginning and at the end of the rising part of the pulse. 
One can see that in Fig. OJa) the analytical curves show 
oscillations as well as the numerical curve; however, the 
amplitude of oscillations is much larger for the analytics 
(P e goes to zero at some values of r, in contrast to the nu- 
merical results). Neglecting oscillations, we see that the 
numerical curve is approximately in the middle between 
the two analytical curves (on the logarithmic scale). In 
contrast to that, in Fig.[3fb) the numerical curve is close 
to the lower analytical curve (and shows some oscilla- 
tions, which are not well pronounced in the analytical 
results), while in Fig. E[c) the numerical curve is close 
to the upper analytical line. As follows from the results 
in Fig. [3 even though Eq. i)28|) can be used for a simple 
estimate of the nonadiabatic error, the accuracy of this 
estimate can easily be worse than one order of magni- 
tude. [The ratio between the upper and lower analytical 
estimates in Figs.EJb) and[SJc) is crudely r£, where r w 
is the ratio of frequencies o>io before and in the middle of 
the pulse; and even though r u ~ 1.6 is not a big number, 



FIG. 5: Dependence of the measurement error P e on the pulse 
duration r calculated either numerically (solid lines) or using 
the simple formula l|28|l . in which the spectral component is 
taken at the frequency Wio corresponding to either beginning 
(dotted lines) or end (dashed lines) of the measurement pulse. 
The pulse is (a) linear, g(6) = 9, (b) of the shape g{0) = 
1 — cos 4 (7r#/2), or (c) of Gaussian shape. For the Gaussian 
pulse r is defined as the r.m.s. width of the pulse. 

is already more than an order of magnitude.] Notice 
that P e for the Gaussian pulse [Fig. Etc)] scales with r 
as t~ 4 , and not in a Gaussian way, as might be naively 
expected since Fourier transform of a Gaussian is a Gaus- 
sian. The reason is that in our case we should use only 
the rising half of the Gaussian pulse, for which Fourier 
transform scales with frequency only as lu~ 2 . 

The numerically calculated dependence P e (r) for sev- 
eral pulse shapes [shown in Fig.^a)] is presented in Fig. 
EJb). We consider the following shapes g{9): 

6, sin(7r0/2), 1 - cos(tt0/2), 

sin 2 (7r6»/2), sin 4 (7r0/2), 1 - cos 4 (7r0/2). (39) 

The long-r behavior of the curves in Fig. Elb) can be 
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FIG. 6: (a) Several shapes of the measurement pulse and 
(b) the corresponding dependences P e (j). The pulse shapes 
are g(6) — 9 (thin solid line), sm(ir8/2) (thin dotted line), 
1 — cos(7r6'/2) (thin dashed line), sin 2 (7r#/2) (thick solid line), 
sin 4 (7r#/2) (thick dashed line), and 1 — cos (n8/2) (thick dot- 
ted line). 



checked to satisfy Eq. I|31l) . In particular, at long r the 
error decays approximately as t^ 2 -? where j [defined right 
after Eq. (|31|l as the lowest order of nonvanishing deriva- 
tives at the pulse endpoints] is equal to 1 or 2 for the 
shapes shown in the first or second lines of Eq. (|39|l . 
respectively (thin and thick lines in Fig. EJ. By design- 
ing pulses to be even smoother at the endpoints (having 
larger j) the decay of P e with increasing r can formally 
be made even faster. However, this does not make much 
practical sense because for long r the error is very small 
anyway, and because the derivation of Eq. (|31[1 required 
an assumption that the pulse shape is described by an an- 
alytical function; this assumption can hardly be applied 
to an experimental situation. 

From experimental point of view the most interesting 
range of P e is 10 -4 — 10 -2 (recall that the real measure- 
ment error can be up to Q e — \fP~l) . We see that in this 
range P e still significantly depends on the shape of the 
pulse. Therefore, for a fast and reliable measurement the 
choice of a proper pulse shape can be really important. 
Notice that the oscillating behavior of P e (T) dependence 
for some shapes can in principle be used to minimize 
nonadiabatic measurement error by choosing specific du- 
ration of the measurement pulse. 



IV. DISCRIMINATION BETWEEN QUBIT 
STATES BY TUNNELING 

In this section we neglect the nonadiabatic error dis- 
cussed above and consider the measurement error arising 
during the tunneling stage of the measurement process. 
Ideally, the system in state |1) should tunnel from the 
left to the right well of the qubit potential (Fig. 1), while 
state |0) should remain in the left well. However, since 
the ratio of tunneling rates Ti/Tq for the two states is 
finite, discrimination between the states by tunneling is 
not complete: for a too short measurement pulse there is 
a chance that the state |1) remains in the left well, while 
for a too long pulse there is a chance that the state |0) 
tunnels out. 

To take tunneling into account we consider the exact 
potential @ and still assume qubit parameters of Eq. 
. Figure shows the dependence of the energy spec- 
trum on the applied flux cf>. The flux is normalized by 
the critical value 4> c at which the barrier between the two 
wells disappears (for our parameters <p c = 5.43). Zero of 
the energy is chosen at the bottom of the left well for 
(j> < 4> c and at the point of inflection of the potential 10 
S c = 7r/2 + arcsin(l/A) for cf) > <f> c . The dashed line shows 
the barrier height, which satisfies quite well the cubic- 
potential approximation2ii2iii AUi cx (1 — 0/</> c ) 3 / 2 . This 
line separates states above the barrier which are delocal- 
ized and states under the barrier which are localized ei- 
ther in the left well or in the right well. In Fig. almost 
horizontal branches are the left-well states. They are de- 
noted on the left side of the frame as |fc) (Ar = 0,1,...), 
where k enumerates the left-well states. The branches 
with a steep slope below the barrier top correspond to 
the right-well states. At <fi > <fi c the potential represents 
a single well. Notice the (anti)crossings between energy 
levels belonging to the left and right wells. In the case 
when the energy relaxation and dephasing are sufficiently 
weak so that the corresponding broadening of the right- 
well levels is smaller than their energy separation, the 
level crossings correspond to enhanced tunneling from 
the left well. 

For simplicity, let us first neglect the level discrete- 
ness in the right well and calculate the tunneling rates 
Ti and Tq for states |1) and |0) using WKB approx- 
imation. Then Tk = f a (k)D(k), where WKB factor 

D(k) = exp[-(2/7i) ^2m(U(S) - E k )dS] is deter- 
mined by the integral between classical turning points 
Sti and 6t2 (we use notation Ek for the energy of the 
qubit state \k) to distinguish it from notation used 
in Sections II and III). For the attempt frequency f a (k) 
we use approximation^ i" a (fc) = (uji/2tt)[(2tt) 1 / 2 /k\)[(k + 
l/2)/e] fe+1 / 2 , where uji is the left-well plasma frequency. 
(Even though this formula has been derived 2 ^ for the case 
when there are many levels in the well, and so it is not 
quite accurate in our case, the inaccuracy is not very sig- 
nificant for us.) Figure |H1 shows the tunneling rates Ti 
and Tq as functions of 4>/4> c - The lines end when the cor- 
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FIG. 7: Structure of the energy levels E n (in GHz) as a func- 
tion of external flux (f> (in units of critical flux <f> c = 5.43). The 
levels with n between 168 and 218 (counted from the right- 
well bottom) are shown. The dashed line shows the barrier 
height AUi. The vectors \k) on the left enumerate the states 
localized in the left well. 




FIG. 8: The WKB tunneling rates Ti and T for the states 
|1) and |0) as functions of the applied flux cf>. 



responding levels reach the barrier top. Small steps near 
the ends of the lines are due to level anticrossings. One 
can see that the ratio Ti/T$ decreases with increasing 
flux cf> and is slightly over 10 2 in the practically interest- 
ing range where Ti > 10 s s . 

Since the applied flux <j) changes during the mea- 
surement pulse, the tunneling should be integrated over 
the pulse duration. For definiteness let us consider the 
parabolic measurement pulse 

<P(t) = 0o + 4(0x - <Po)(t/r)(l - t/r) (0 < t < r), (40) 

so that the flux (j)(t) increases from 0o to the maximal 
value 4>\ at t = r/2 and then decreases back to 0o a * 
t = t. (Notice a different definition of the pulse du- 
ration r compared to the definition used in Sec. III. in 
which we considered only half-pulses.) Figure [5] shows 
the probabilities of tunneling during the measurement 
pulse (switching to the right well) P Sj i and P Sj o for levels 
|1) and |0), 
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FIG. 9: The probabilities P s .i (thick lines) and P 3t o (thin 
lines) of tunneling to the right well during the measurement 
pulse, starting from the states |1) or |0), correspondingly, as 
functions of the maximum flux 4>\ during the pulse. Calcu- 
lations are based on the WKB approximation. Initial flux is 
4>o/4>c = 0.94 (N = 5). The pulse duration is r = 2 ns (solid 
lines), 10 ns (dashed lines), or 50 ns (dotted lines). 



as functions of the maximum flux 4>\ (normalized by crit- 
ical flux (j) c = 5.43); initial flux O = 5.09 (<V0 C = 0-94) 
corresponds to N — 5; the pulse durations are t = 2 ns, 
10 ns, and 50 ns. Notice that Tk is not well-defined when 
the state goes over the barrier, while we need some value 
of Tk for integration in Eq. 141|) ; in this case we still use 
the definition Tk = f a (k)D(k) but assume D(k) = 1. 
The curves in Fig. remain at practically zero level 
for small pulse amplitudes and saturate at 100% level 
for large enough pulse amplitude (the "S-curve" shape). 
(The switching probability P S; o for t = 2 ns does not 
fully approach 100% because the pulse duration is too 
short.) The most important observation is that the flux 
shift between the curves for P s i and P Sj o (for the same 
t) is sufficiently large to reliably distinguish states |1) 
and |0). However, the measurement error, which can be 
defined as (P Sl o + 1 — Ps,i)/2, is finite for any <j>i, and 
in the optimal point is on the order of 10~ 2 . (Actu- 
ally, our method is not quite accurate for the top part 
of the S-curves because of significant level anticrossings, 
which leads to visible kinks on the curves in Fig. [5] and 
inaccuracy of the measurement error calculation.) One 
can see that the minimal measurement error improves 
(decreases) with increase of the pulse duration t; even 
though the separation of the two S-curves decreases with 
increase of r, they become sharper, leading to better dis- 
crimination between states |1) and |0). 

The WKB approximation does not take into account 
the level discreteness in the right well. To take it into ac- 
count, we have performed the full quantum-mechanical 
simulation of the evolution Q) (JSJ) during the measure- 
ment pulse i|40|) ■ Since the level broadening is determined 
by energy relaxation (and dephasing), which is very dif- 
ficult to simulate exactly, we have used a simpler model. 
We assume that the levels in the right well and above 
the barrier have a finite lifetime (complex energy) , while 
levels in the left well do not decay. The decay rate for 
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the nth level counting from the bottom of the right well 
is chosen as 7 = n/T\, where T± is the energy relaxation 
time; the corresponding width of the broadened level is 
hj. Since there are many levels in the right well and we 
are interested only in levels close to the barrier top, the 
level index n can be replaced with the number of levels 
N r in the right well; for the parameters J7J which we use, 
N r = 174 at TV = 5. [Actually, there are some additional 
details in the algorithm. For example, to prevent unphys- 
ical decay of the left-well state during the level crossing, 
we smoothly (parabolically) suppress 7 when the energy 
distance to the nearest level becomes less than 1/5 of the 
normal level separation in the right well. The numerical 
criterion for a left-well level is that its energy is in be- 
tween the left well bottom and the barrier top, and also 
the maximum of the wavefunction is to the left of the 
barrier top.] 

Solid lines in Fig. show the numerically calculated 
probability P s k of switching to the right well during the 
measurement pulse (starting with N = 5, 4>/4> c = 0.94) 
in the case when initially either the left-well ground state 
(k = 0) or the first excited state (k = 1) is occupied. As- 
sumed pulse duration is r = 2 ns, while T\ — 60 ns. 
Switching probability is defined as P s = 1 — Pi, where 
Pi is the total population of the left-well levels after 
the pulse. One can see that the numerical dependence 
-P«,fc(0) is generally similar to the corresponding WKB 
result (dashed lines); however, the numerical curves show 
significant oscillations due to level discreteness in the 
right well. Notice that the amplitude of oscillations in 
this case is mainly determined not by the width of the 
levels (which is still relatively small in this example), 
but by continuous sweeping through the comb of levels 
during the pulse. (For our parameters the level broad- 
ening becomes larger than the level spacing at T± < 5 
ns.) Similar oscillations in the dependence P s ,k(4>i) have 
been observed experimentally^ though for a qubit with 
a smaller number of levels in the right well (N r ~ 30, in 
contrast to N r ~ 170 in the example shown in Fig. imp . 
that leads to a smaller level broadening ^7 ~ HN r /Ti 
and therefore more pronounced effects of the right-well 
level discreteness. 

Figures ISlllOl show that the quality of tunneling dis- 
crimination between states |1) and |0) depends signifi- 
cantly on the choice of the measurement pulse amplitude 
(maximum flux cf>i), and at the optimal point the mea- 
surement error (P s ,o + 1 — Ps.i)/2 is on the order of 10~ 2 . 
This error decreases for longer measurement pulses and 
can be further reduced by changing the qubit parameters. 



V. LEFT WELL REPOPULATION 

A steady progress in fabrication of phase qubits has re- 
sulted in a significant decrease of the dissipation ratei2£ 
This trend will obviously continue because of the require- 
ment of very low decoherence in quantum computation. 
Therefore it is of interest to study the measurement pro- 
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FIG. 10: The switching probabilities P s ,o (thin lines) and 
P Si i (thick lines) as functions of the maximum flux <f>i, cal- 
culated numerically (solid lines) and using WKB approxima- 
tion (dashed lines). The measurement pulse is of parabolic 
shape with duration r = 2 ns, starting and ending at flux 
4>o/4>c = 0.94 (N — 5). Oscillations of the solid lines re- 
flect level discreteness in the right well. (At N — 5 there are 
N r = 174 levels in the right well.) 



cess for a phase qubit with a very long energy relaxation 
time T\. Let us we consider the limit of negligible dissi- 
pation, T\ — 00. As we will see below, in this case the re- 
sults are qualitatively different from the results discussed 
in the previous section. The most important effect is that 
the switching probability P S: k does not approach 100% 
even for measurement pulses with maximum flux <\>\ well 
above the critical value 4> c , in contrast to what is shown 
in Figs. 131 and ITU1 

To demonstrate this effect, let us assume that initially 
the system is in the left- well ground state |0) and simu- 
late the full quantum evolution Q-0 during the mea- 
surement pulse (|40|l . Figure 1111 (the main frame) shows 
the populations P n (t) of levels n (levels are counted from 
the right-well ground state) with 169 < n < 204 for the 
measurement pulse with initial flux <po = 0.94(/> c (corre- 
sponding to iV = 5), maximum flux <pi = l.01(f> c , and 
total duration r = 2 ns. The inset shows the structure 
of energy levels E n (for 168 < n < 218) as a function of 
time during the first (increasing) half of the pulse (the 
only difference between Fig. [7J and the inset of Fig. ^] 
is the converted horizontal axis). The state |0) initially 
corresponds to n = 169. 

The population of the qubit state |0) can be seen as 
the upper envelope dependence of P n (t) for t < 0.63 ns; 
similarly, lower envelope curve at t < 0.4 ns corresponds 
to the state |1). At t < 0.4 ns the state |0) is sufficiently 
deep in the left well, so that the level splittings in the 
avoided crossings are too small to cause transitions from 
this level to the right well (notice though the change of 
the level number corresponding to |0) after each level 
crossing). As a result, at t < 0.4 ns the populations of 
levels |0) and |1) vary only due to the nonadiabatic effects 
(see Secini. 

At 0.4 < t < 0.63 ns there is an appreciable probability 
of a transition to the right well at each crossing of the 
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FIG. 11: The level populations P n (169 < n < 204) as functions of time t for the measurement pulse with duration r = 2 ns, 
starting and ending at <f>o/4> c = 0.94 (N = 5) with maximum flux in the middle <j>i/(j>c = 1.01. The curves for 169 < n < 197 
are thin solid lines, the curves for 198 < n < 204 are denoted by the value of n. Inset: Energies E n (168 < n < 218) versus 
time t for the rising part of the pulse (similar to Fig. 0. 



qubit level |0) with a right-well level. In this time inter- 
val one can still interpret the envelope of the populations 
of the crossing states as the population of state |0), de- 
spite the level |0) becomes above the barrier at t > 0.5 
ns (see the inset). We have checked numerically that the 
transition probabilities at the level crossings agree with 
the Landau-Zener formula 18 . At t < 0.5 ns, when the 
state |0) is still below the barrier, the transitions to the 
right well are due to tunneling. At t > 0.5 ns the level 
splittings increase and become comparable to the sep- 
arations between adjacent right-well levels; this means 
that generally more than two levels are coupled simul- 
taneously. At t > 0.63 ns the left well disappears (see 
the inset), and the initial population of level |0) becomes 
distributed between several states in the resulting single 
well. 

During the time interval when the potential @ has a 
single-well shape (0.63 ns < t < 1.37 ns), one can see a 
rather significant and complicated redistribution of the 
level population. Notice that the most populated levels 
in the single well are those with 198 < n < 204, i.e., with 
the quantum number significantly higher than the initial 
n = 169. The solid line in Fig. 1121 shows the population 
distribution at t = 1 ns (middle of the pulse); one can see 
that the population is noticeable only for n > 192, which 
are the states delocalized (above-the-barrier) before and 
after the pulse. 

After t = 1 ns the pulse starts to decrease, and the 
potential @ starts to return to its initial shape. The 
levels localized in the left well start to appear again and 
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FIG. 12: The level populations P n at t = 1 ns and t — 2 ns 
for the process shown in Fig. 1111 (cbi I6r = 1.01, t — 2 ns). 
The data are shown by dots and squares, connected by lines 
as guides for the eye. 



cross with the populated right-well levels which are mov- 
ing up. As a result, the qubit state |0) partially recovers 
its population (in Fig. 1111 it is seen as an approximately 
horizontal line after 1.6 ns), and also a number of higher 
levels in the left well become populated. As the dashed 
line in Fig.li"2lshows. after the end of the pulse an appre- 
ciable population is acquired by states |0), |1), and |2) of 
the left well and also by a number of delocalized states 
with n > 190. 

In short, during the first half of the pulse, right-well 
and delocalized levels go down in energy and cross with 
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left-well levels, thus gaining population, while during the 
second half of the pulse, right- well levels move up and 
cross again with left-well levels, resulting in a partial re- 
population of the initial state and an excitation of higher 
levels. Naively, one could expect that after such strong 
measurement pulse the population should be transferred 
from the left well to the right well. However, we see that 
without dissipation the population actually goes mainly 
to highly excited delocalized states. Assuming a very 
small but finite dissipation after the pulse, we conclude 
that a populated delocalized state will eventually end up 
either in the left or in the right well. Therefore, there are 
two contributions to the left well repopulation: direct re- 
population (due to level crossings) and relaxation from 
delocalized states. 

After the pulse the probability for the system to be to 
the right of the barrier top is 



P r {t) 



\V(6,t)\ 2 d6 



(42) 



(where S m is the barrier-top position), and generally de- 
pends on time due to interference of occupied delocalized 
levels. We define the probability of escape (switching) to 
the right well as the time average of P r (t). Thus, the 
switching probability is given by 



P, 



\M$)\ 2 dS, 



(43) 



where P n is the population of level n after the pulse, and 
the integral is the probability for the system in state n 
to be in the right well. This integral equals zero (one) 
for states localized in the left (right) well and is between 
and 1 for delocalized states. 

The solid line in Fig. I13f a~) shows the dependence of 
the switching probability P s .o (starting from state |0)) 
on the maximum flux <pi for r = 2 ns and </>o/0 c = 0.94. 
The switching probability practically vanishes at <fi\/ <j) c < 
0.98 when the initial state |0) remains sufficiently deep 
inside the left well and cannot tunnel (similar to the re- 
sults shown in Figs. I§1 and 110(1 . However, in contrast to 
Figs. I$l and 1101 an increase of the pulse amplitude leads 
to significant and irregular oscillations of P S! i(^>i). Most 
importantly, the switching probability remains to be sig- 
nificantly less than 100% even for 4>i exceeding C . 

The solid line in Fig. I13f b) shows the switching prob- 
ability P Sj o as a function of the pulse duration r for 
(j>i = 4> c - The dependence is even more irregular than in 
Fig. I13f ah The irregular oscillations seen in Figs. I13f ai 
and !13f bi can be related to the complicated time depen- 
dence of the level populations seen in Fig. ITTI 

Besides calculations for the qubit parameters shown in 
Eq. Q), we have also considered a different set of qubit 
parameters: 



C = 790 fF, L = 0.720 nH, I = 0.764 /xA. 



(44) 



The main difference between the two sets is a significantly 
different number of levels in the right well. For the initial 
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FIG. 13: The switching probability P s ,o for the initial state 
|0) [see Eq. I)43|l ] (a) versus the maximum flux <f>\ (normalized 
by the critical flux (f> c ) for r = 2 ns and (b) versus the pulse 
duration r for <f>i — cj> c . Solid lines are for our usual set 
Q of qubit parameters; dashed lines are for the qubit with 
parameters 144H . having much smaller number N r of discrete 
levels in the right well. 



condition N = 5, the number of levels N r in the right 
well is N r — 174 for our main set of parameters J7J and 
N r = 30 for the parameters (|44() . The eigenenergies of 
the Hamiltonian as functions of flux <fi for the parameter 
set <|44[) are shown in Fig. ^] (flux is normalized by the 
critical flux (j) c — 3.55). 

In spite of a significant difference of the energy level 
structure, the behavior of the switching probability P s q 
shown in Figs. I13f a) and I13f b) by dashed lines for the 
parameter set (|44f) is qualitatively similar to the previ- 
ously discussed behavior for the set J7J) (solid lines) . No- 
tice, however, that in a real experiment the effective time 
7 _1 ~ Ti/N r spent on a right-well level before relaxation 
is significantly larger for the set l|44l) and therefore the 
discussed here effects due to insufficiently fast relaxation 
will be more pronounced. 

It is important to notice that while the large irregu- 
lar oscillations of switching probability P s in Fig. ^] is 
a result of quantum interference and therefore will be 
significantly suppressed by a weak relaxation, the fact 
that P s does not approach 100% even for 4>i > <j) c [as for 
curves in Fig. I13f a)] is a much more robust, essentially 
classical effect. Our results show that a very significant 
repopulation of the left well may occur if the dissipation 
is insufficiently fast. In other words, sufficiently strong 
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FIG. 14: Energies E n (in GHz) as a function of flux (f> (in 
units of 4>c = 3.55) for the qubit with parameters 144H corre- 
sponding to N r = 30. The levels with 24 < n < 57 are shown. 
The dashed line corresponds to the barrier height AUi. 

energy dissipation is a necessary requirement for a good 
measurement fidelity. Conversely, the repopulation effect 
imposes a lower limit on the pulse duration for a given 
dissipation rate. 

Let us estimate the condition necessary to eliminate 
the left-well repopulation. For simplicity, we assume the 
pulse to be almost rectangular, with the rise and fall 
times being significantly shorter than the time r m , dur- 
ing which the pulse is at its maximum. Let us denote by 
Ui the level number (counting from the right-well ground 
state), corresponding to the right-well level close to the 
top of the barrier for maximum flux bias. Assume that 
this level is populated immediately after the pulse rise. 
Then during r m the energy relaxation will lead to the 
decrease of n-number corresponding to the mostly pop- 
ulated level. Neglecting the "energy diffusion" of the 
populated levels, we estimate the n- number at time r m 
as v! i ~ rii — 7T m , where 7 is the transition rate between 
adjacent levels. We neglect the weak dependence of 7 
on the level number. Let rif be the first level above the 
barrier top after the pulse is over. If we neglect a change 
of •n! i during the falling part of the pulse, the condition 
that there is no return to the left well is •n! i < rif, which 
means r m > (rii — U/)/7- Using an estimate 7 ~ rii/Ti, 
we obtain 



This condition is actually a lower bound; the required 
pulse duration r TO can be even longer because of not im- 
mediate population of the level rii, widening of the range 
of populated levels in the course of relaxation ( "energy 
diffusion"), and possibility of repopulation due to tun- 
neling. For the case shown in Figs. ^] an( i IS! we nn d 
m = 205, rif = 181, yielding r m > 0.13 T x . 

The necessity of sufficiently long measurement pulses 
thus seems to be an important requirement for the design 
of quantum gates based on phase qubits. It also makes 
harder the solution of the cross-talk problemSii£ for si- 
multaneous measurement of several qubits, which argues 



in favor of designing adjustable coupling between qubits. 



VI. CONCLUSION 

In this paper we have studied the behavior of a flux- 
biased phase qubit in the process of its measurement (us- 
ing analytical approaches and numerical solution of the 
Schrodinger equation) and analyzed several mechanisms 
leading to measurement errors. 

First, we have studied nonadiabatic errors (Sec. Ill), 
which occur during the rising part of the measurement 
pulse (before the tunneling stage) due to hnite duration 
of the pulse, leading to the transitions between qubit 
states |1) and |0). We have developed simplified analyt- 
ical approaches with several levels of accuracy [see Eqs. 
(1251 . 1271) . and Q28p] and compared them with the numer- 
ical results (Sec. IIII B|) . The nonadiabatic error generally 
decreases with increase of the pulse duration r; however, 
this dependence may exhibit significant oscillations. The 
numerical value of the error depends significantly on the 
pulse shape (Figs. |5] and |BJ), thus showing importance 
of a proper pulse shape design. For pulse duration over 
few nanoseconds the nonadiabatic error is typically suf- 
ficiently small to be considered negligible. 

Another type of measurement error is due to incom- 
plete discrimination between states |1) and |0) during 
the tunneling stage of the measurement process (Sec. IV). 
For typical qubit parameters^ the ratio of WKB tunnel- 
ing rates T\/Tq for states |1) and |0) is few times 10 2 and 
decreases with the applied flux bias (i.e. when faster mea- 
surement is required). Therefore, the minimized (over 
the maximum applied flux during measurement) mea- 
surement error is smaller for longer pulses (Fig. . The 
level discreteness in the right well of the qubit potential 
[Fig. [Hb)] leads to oscillations in the dependence of the 
switching probability on the measurement pulse ampli- 
tude (Fig. llOfl . A typical value of the minimized error 
due to incomplete discrimination between the states by 
tunneling is on the order of 10 -2 for typical qubit pa- 
rameters of present-day experiments. 6 Change of qubit 
design may lead to further reduction of this type of mea- 
surement error. 

In Sec. V we have analyzed the quantum evolution of 
the qubit during its measurement in the case of complete 
absence of energy relaxation, and found that the energy 
relaxation plays a very important role in the process of 
qubit measurement. In the case of insufficiently fast dis- 
sipation (characterized by energy relaxation time T{) the 
measurement error can be caused by the repopulation 
of the left well of the qubit potential [Fig. db)] after 
the tunneling stage of the measurement process. A sim- 
ple estimate l|45[) for the repopulation effect shows that 
reliable measurement requires relatively long measure- 
ment pulses; in the analyzed example the pulse should 
be longer than ~ 0.1 T\. This result may be quite impor- 
tant for the design of the quantum gates based on phase 
qubits. Some alleviation of the problem may be achieved 
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by using a different shape of the measurement pulse, so The work was supported by NSA and DTO under ARO 
that the barrier height of the qubit potential after the grant W911NF-04- 1-0204. 
pulse is significantly lower than before the measurement. 
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